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Abstract 

A gauge-invariant, linear cosmological perturbation theory of an almost homogeneous and 
isotropic universe with dynamically evolving Newton constant G and cosmological constant A 
is presented. The equations governing the evolution of the comoving fractional spatial gradients 
of the matter density, G and A are thus obtained. Explicit solutions are discussed in cosmologies, 
featuring an accelerated expansion, where both G and A vary according to renormalization group 
equations in the vicinity of an ultraviolet fixed point. Finally, a similar analysis is carried out in 
the late universe regime described by the part of the renormalization group trajectory close to the 
gaussian fixed point. 
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I. INTRODUCTION 



The subject of cosmological perturbations continues to attract much attention because 
it is an essential step in understanding any theory of cosmological structure formation. The 
simple idea that the observed structure in the universe has resulted from the gravitational 
amplification of small primordial fluctuations works remarkably well and it must be discussed 
in any theory of relativistic cosmology. 

In inflationary cosmology the presently observed structures in the universe are generated 
by quantum fluctuations during an early de Sitter phase. The subsequent evolution is 
classical and depends on the interplay between pressure forces, the rate of growth of the 
expansion factor and on the content and the nature of the, yet unknown, dark matter. 

Quite recently, an alternative scenario was proposed. It is based on the assumption that 
the gravitational couplings are subject to renormalization group (RG) flow, and reach a fixed 
point when the energy goes to infinity. Evidence for this has been found first in the Einstein- 
Hilbert truncation [1-7], followed by extensions to four-derivative gravity [8-10] and to f(R) 
gravity [11-14]. Further extension have been discussed in the conformal truncation, where 
one freezes the spin two degrees of freedom [15, 16], and in bimetric truncations [17]. 

In [18] the RG scenario was applied to a Friedmann-Robertson- Walker cosmological 
model, assuming that the relevant energy scale upon which the couplings depend is the 
Hubble parameter H(t). This gave rise to a "RG improvement" of Einstein equations that 
was found to be responsible for a phase of accelerated expansion of the scale factor. RG 
scenario was thus able to trigger a mechanism of "inflation" without the introduction of an 
ad hoc scalar field. 

In [19] the evolution of gauge- invariant perturbations, according to the formalism pio- 
neered by Hawking [20], further extended by Ellis and Bruni [21], Jackson [22] and Zim- 
dahl [23], was studied in a homogeneous and isotropic background universe, deep inside the 
late universe regime (where the authors postulated the existence of an infrared fixed point). 
It was found that all the relevant perturbations decay with a power law of the cosmic time or 
stay constant, i.e. the IR fixed point cosmology is stable under small deviations from perfect 
homogeneity and isotropy. The aim of the present paper is to investigate the evolution of 
cosmological perturbations in the vicinity of the UV fixed point of the RG trajectory, close 
to the above mentioned accelerated phase of the universe. 
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It is important to point out the fact that, in [19], the stress-energy tensor of the perfect 
fluid filling the universe was assumed to be covariantly conserved, an hypothesis that only 
works well for low values of the cutoff scale, when the entropy exchange between the matter 
sector and the coupling constants is suppressed (see [18] for details). Moving backwards 
along the RG trajectory up to the Planck mass scale, the suppression breaks down and the 
sudden variation of the cosmological constant A generates a flow of energy from the geometry 
into the fluid. As a consequence, the conservation of the stress-energy tensor gets spoiled 
and one can only count on the more general contracted Bianchi identity of the Einstein 
tensor. 

The paper is structured as follows: in section II the main features of the RG improved 
background cosmology are briefly reported, with a closer look at the very early stages of its 
evolution. In section III the covariant approach to the evolution of perturbations is described 
and is applied to the RG-driven inflationary background, getting power law solutions that 
will be discussed in section IV. Finally, in section V the same analysis is applied to the 
subsequent phase of quasi-classical evolution predicted by the RG improved gravitational 
theory. 

II. RENORMALIZATION GROUP DERIVED BACKGROUND 

The idea of using the RG flow equation in gravity is borrowed from particle physics 
where the RG improvement is a standard device in order to add, for instance, the dominant 
quantum corrections to the Born approximation of a scattering cross section. However, 
instead of improving solutions, in [24, 25] the more powerful improvement of the basic 
equations has been discussed. The main results for an homogeneous and isotropic universe 
will now be briefly reviewed. 

It will be assumed that there is a fundamental scale dependence of Einstein-Hilbert 
coupling constants which is governed by an exact RG equation for a Wilsonian effective 
action whose precise nature needs not to be specified here. Regardless of the technical details 
of the chosen renormalization scheme, there is a rather general fact about renormalization 
group improved actions that is worth to be stressed. A renormalized action is composed 
by a sum of operators that depend both on the fields, that we generically call 0, and their 
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momenta p, coupled to constants depending on a renormalization scale k 

r fc = X>(*0<W,p)- (!) 

i 

As such scale is typically identified as k ~ p, one gets that 

r^Yldiipm&p) = Yl°i^p) (2) 

i i 

where the new operators get additional dependences on the momenta, that are rather generic 
and can even include nonlocality. This means that the renormalization group itself is a kind 
of shortcut to take into account more terms than the ones actually present in the action. A 
more exhaustive discussion about this topic can be found in [26]. 

Getting back to the Einstein-Hilbert framework, at a typical length scale £ or mass scale 
k = l~ x the coupling constants assume the values G(k) and A(k), respectively. In trying 
to RG-improve Einsteins equation the crucial step is the identification of the scale which is 
relevant in the situation under consideration. In cosmology, the postulate of homogeneity 
and isotropy implies that k can depend on the cosmological time only so that the scale 
dependence is turned into a time dependence: 

G = G(k(t)) , A = A(k(t)) . (3) 

In [18, 27] this dependence was implemented by setting such energy scale to be propor- 
tional to some effective measure of the curvature of the universe, i.e. the Hubble rate: 

k(t)=tH(t) = t- (4) 
a 

where £ is a 0(1) parameter that is kept free for the time being. For the sake of completeness, 
it has to be pointed out that other choices are possible, as for exemple the cutoff k ~ 1/t 
proposed in [24, 26, 28] or the k ~ l/a(t) in [29, 30], while in [31-33] the cutoff dependence 
was adjusted to fit the energy-momentum conservation of ordinary matter. 

Later on, when dealing with perturbations around the Friedmann-Robertson- Walker 
(FRW) background, one will have to give up perfect homogeneity and isotropy, so equa- 
tion (3) will generalize to 

G = G(k(x»)) , A = A(k(x»)) . (5) 

In this context H is not defined anymore, and in section III it will be identified k = £9/3, 
where 9 is the expansion rate of a congruence of world lines of the cosmological fluid. In the 
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FRW case 6 = 3H, so there is consistency with the choice made in (4). As said before, the 
cutoff scale is meant here to be identified with some measure of the local curvature, and gets 
a coordinate dependence only because the curvature happens to be coordinate-dependent. 
From this point of view, there is no conceptual difference between the choice k = k(H(t)) 
in a FRW universe and the one k = k(6(x^)) in a quasi-FRW one. 

For what regards the homogeneous and isotropic background, in the following the main 
results of [18] are briefly resumed. With G(k) = g{k)k~ 2 and A(k) = X(k)k 2 , g(k) and X(k) 
being the dimensionless couplings, one has 

G(t) = , A(t) = A(£ H(tm 2 H(t) 2 • (6) 

Then, provided H(t) ^ 0, the cosmological evolution equations for a fluid with equation of 
state p = wp can be cast in the following form: 



H(t) = ~(l+w)H(t) 2 
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(7) 
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tit) 



H{tf . (8) 



8ng(ZH(t)) [ 3 

In [18] it is stressed that for k — > oo the RG flow spirals onto a Non-Gaussian Fixed Point 
(NGFP), whose values of the dimensionless coupling constants are denoted by g* and A*. In 
that paper the spiral is approximated to a constant-^ and A trajectory, so the differential 
equation (7) reads 

H = -a^H 2 (9) 

with the constant 

2 _ 2 

a = 3(1 + w)(l - Ke/3) - 3(1 + w)(l- Sl* A ) (10) 

Q* A being the NGFP value of n A (£H(t)) = £ 2 A(£ H(t))/3. 

Fixing the constant of integration such that this singularity occurs at t — 0, the unique 
solution to (9) reads 

H(t) = j (11) 

which integrates to a(t) oc t a . Using il* A as the free parameter which distinguishes different 
solutions, the fixed point cosmologies are characterized by the following power laws: 



a(t) = a t 
5 



tit) 



t 

9ng,K(l+w) A (l-n* A ) 3 

3(?,A,(1+«;) 2 (1-^) 2 + 
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G(t) 



(12) 



A(t) 




Note that the RG data enter the solution (12) only via the universal product g*A* [4, 5]. By 
the analysis of the deceleration parameter 



it can be seen that, in a radiation-dominated regime, the universe undergoes an accelerated 
expansion if Q* A > 1/2. In section IV we compute gauge-invariant perturbations in a NGFP 
background, which corresponds to the inflationary phase of the RG improved cosmology. 

There is one last thing about the evolution of the background universe that needs to be 
defined: after having left the NGFP, the trajectory flows towards the Gaussian Fixed Point 
(GFP), passing very close to it. This means that the couplings become small, and then the 
flow can be linearized and the /9-functions assume the form 



q = - - 1 = ^[1 + 3w - 3(1 + w)Q* A ] 



a 2 



(13) 
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and the corresponding trajectory can be parametrized as 




(15) 



where At = A (fey) and kr is the value of the cutoff scale for which 



(3 x (k T ) = k—X(k) 



= 



(16) 



k=kT 



and that is called turning point. 

During this epoch, the dimensionful coupling constants behave like 



G(k) = G 

A(jfe) = A + uGk 4 
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(17) 



where v is a constant of order unity whose precise value depends on the renormalization 
scheme, and at very late times they approach their present observed values G, A . This is 
the reason why such regime is denominated quasi-classical. These formulas will prove useful 
in section V, when dealing with evolution of perturbations in the late universe. 



III. GAUGE-INVARIANT PERTURBATION THEORY 

In this section a generalization of the formalism of refs. [20-23] is presented, by allowing 
G = G(x^) and A = A(x M ) to be scalar functions on spacetime. 

A fundamental observer describing the cosmological fluid flow lines has 4-velocity 

vJ 1 = dx^/dr, = -1 (18) 

where r is the proper time along the flow lines. The projection tensor onto the tangent 
3-space orthogonal to is 

hfu, = g^u + UnUv (19) 

with h^Jf v = h^ u and h^ u u u = 0. The covariant derivative of is 

u^u = + a + - ii^u u (20) 

where <j^ u = h a hP v U( a .^ — 9h^ u /3 is the shear tensor, = h a ^ u U[ a; ^ is the vorticity 
tensor, 9 = u^.^ is the expansion scalar and ii^ = u fM . v u u is the acceleration four-vector 
(square and round brackets denote anti-symmetrization and symmetrization, respectively). 
The Riemann tensor is defined by w M;pcr — u n;°p = -^po-^^A and the Ricci tensor R^ u = RJ ' . 
With these conventions, Einstein s equations are 

G^ u = R^ u - ^Rg^u = SttGT^ - Ag^ u (21) 

where A = A(x^) is the coordinate-dependent cosmological constant and G = G(x 11 ) the 
coordinate-dependent Newton constant. As already said, the conservation constraint for 
energy-momentum tensor = is going to be dropped in favour of the less restrictive 
hypothesis 

G^. v = {8nGT^ - kg^). v = . (22) 

A perfect fluid stress-energy tensor has the form T^ v = pu^u" +ph> ML ', so eqn. (22) leads 
to a modified energy conservation 

p + 9(p + p) + ^p+^ = (23) 
7 



and a modified equation of motion 

" + JT- p {'» + lT>-£b) = (24) 

by projecting along and onto the hyperplane orthogonal to u^. 

By a redefinition of the variables it is possible to rearrange the previous equations in a 
shape that is formally identical to the usual stress-energy conservation: defining the quan- 
tities 

n=G P+ l, V.Gp-l (25) 
equations (23) and (24) get recast in 

n + 6(1l + V) = (26) 

and 

u» + = (27) 

with the new variables playing the role of some effective p and p. 

Combining the Einstein equations and the (20) one gets the well known Raychaudhuri 
equation, that in the language of 1Z and V is written as 

6 + \d 2 + 2(a 2 - u 2 ) - + Au(K + 3V) = (28) 

where 2a 2 = a^a^ and 2w 2 = u^u^. The term -u M ;M can be rewritten as [23] 

^ = - kV (^) + k ^^P- (29 » 

the scalar /C is defined as 

K = 2a 2 - -9 2 + 16tt^ (30) 

and it is possible to show that for zero vorticity (u^ u = 0) it coincides with the Ricci scalar 
of the 3-dimensional hyperplane everywhere orthogonal to . An auxiliary length scale 
S(t) is introduced as the solution of the equation 

H- 

Suitable quantities useful to characterize the spatial inhomogeneities of density, pressure 

and expansion should be, respectively, 

Sh\p. u Sh\p. v 
p+p p+p * 
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and to characterize spatial gradients of G and A can be introduced the dimensionless quan- 
tities 

as it is already done in [19]. However, as in the present case we are working with some 
adjusted density and pressure, the most useful quantities turn out to be 

Sh\K. u Sh\V- v 

One can now write, using (20), (23) and (24), 

h\(Sh\n ;X y = ^Qh\n, v + sh\v, v + sh\{n., v y = 

= -S6h\1l, v - Sh\6, v (K + V)- Sh\n ;u (a\ + u\) (35) 

where in the second line the fact was used that (TZ-^)' = — u v TZ- V . 

The expression (35) can be rewritten as a differential equation involving the adjusted 
gradients (34) 

h\V v + ^pV, + + a%)V v + t, = (36) 
and following a similar procedure one can also obtain an equation for 

h\i v + en, + ^et, + {u\ + a\)t u + 2Sh^(a 2 - u 2 ),,+ 

-Sh\(u v . v ). x + 4n(n + V){V, + 3i7„) = (37) 

where Raychaudhuri equation (28) was used to extract the final expression. 

Notice that both (36) and (37) represent the most general extension of the usual gauge- 
invariant perturbation theory equations in the case of non-constant G and A, in which no 
assumptions are made on their functional form. 

From now on we will focus on a homogeneous and isotropic background universe, for 
which a, v = = ii, = 0, S = a, = 6K/a 2 and 6 = 3H. In such universe the relevant 
equations for the background evolution read 

K. = 2(- V + 8n-Jl) = (3) i? (38) 
3 

9 + 12n(n + V) = {3) R/2 . (39) 

We now consider small perturbations of the motion of the fluid and consequently, up to 
first order in the inhomogeneities, the factors in the equations multiplying the quantities V,, 



11^ and tfj, refer to the background. Under these hypothesis (36) and (37) get the following 
form: 

h\i> v + K^pV, + t, = (40) 

3K V 2 
h\i v + —n^ + 2Ht il + 4Tr(1l + V)V^ + — i7 M = (41) 

where the fact that [23], up to linear order, Sh x (u". v ).^ c± — V 2 i7 M /a 2 was used (V 2 is the 
laplacian on the 3-hypersurface). 

For the sake of semplicity, from now on we will focus on a spatially flat FRW, setting 
K — 0. One also would like to assign to the perturbed quantities an equation of state 
p(p) = wp but, while this translates directly in a relation = wD^, it doesn't happen the 
same for and i7 M . A possible solution is to rewrite the latter quantities in terms of the 
formers 

V» = D„ + rP+-^-^- (42) 
p + p 8n&p + p 

n» = P» + TP--±-^ (43) 

P + P &7T&P + P 

and then apply the usual equation of state. 

There is one last missing ingredient to the recipe: in order to end up with a consistent 
set of equations for spatial perturbations, it is necessary to eliminate the variables V ^ and 
A M in favour of the remaining ones. This is achieved by making explicit use of the cutoff 
identification. It is assumed that the cutoff k appearing in equation (5) is k = £6/3 (which 
reduces to k = £ H in the FRW background). Then, using equations (32) and (33), as long 
as the trajectory keeps close enough to the NGFP, where g ~ g* and A ~ A*, one can write 



Sh\G. v Sh\((6/?,)- 2 g). v 2 2 

r = g = M ' ' Zlf. ~ --Sh u = —t 44) 

" G (6/3)-*g ~ 6 ^ 3H^ 



Sh\A. u Sh u M6/3) 2 X). u 2 2 

where in the last step one could get back to H, at least up to first order, because it is 
multiplied to the perturbed quantity t M . To better understand the meaning of what just 
done it is worth stressing the fact that, by including and A M inside the equations, one 
was treating spatial fluctuations of the coupling constants as some sort of matter fields, 
whose dynamics were governed by the Einstein equations. But the dynamics of G and A 
are described by the RG equations and, by means of the cutoff identification, depend on 
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the expansion rate 6. This means the spatial gradients of the couplings turn out to the 
proportional to the gradient of the expansion rate. 

It is now time to plug all the ingredients inside the system, ending up with the two 
differential equations 

A - 8nGp ■ 8nw(Gp + Gp) - A 



h\D v + 12 T r G(l + w)pH" M 



I 



8%G(l+w)p 11 
HA(8tv P G + A)+ 



6(47rG(l + w)pH) 2 
+32n 2 G 2 p ( (1 + w)p(3(l + w)H 2 + 2H) - 2wHp ) + 



-8ttG ((1 + w)KpH + H(87rwp 2 G - (2 + w)pA + Ap) 



t„ = 



(46) 



+ ( 4*G(1 + W)P + + 12 ^ (1 1 +w)pJJ 



3(1 +w)if 2 -4vrG(l + w)p 



-5M- 



4ttG(1 + w)p- 



(47) 



and, having noticed that /i M A (h v D„y = h u x D u — ii\u v D v = h u x D u , one can combine them 
to form a single, second order differential equation in that for the time being will be 
schematically indicated as 



Ait) h\D v + Bit) h\D v + C(t) = 



(48) 



IV. FIXED POINT EVOLUTION OF PERTURBATIONS 

In the following the general formalism developed above is applied to the NGFP back- 
ground cosmology described in section II. First one needs to expand the perturbed quantities 
in spherical harmonics, that is 

= <j) n (t)^ (49) 

where is the eigenfunction of the spatial Laplacian operator V 2 with eigenvalues — e 2 . 
For every n then equation (48) becomes 

A(t) n (t) + B(t) <p n (t) + C(t) (fit) = (50) 
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where, imposing the background (12) 

A(t \ = 3(1 + w) 

1 ) ~ Qw(i - n* A ) + 2(1 + n* A ) - 3e^(i + w )( 2 n* A - 1) { } 



B(t) 



t 



-1 



C(t) 



(Q* A - 1) (6w(n* A - 1) - 2(^ + 1) + 3e^(l + W )( 2 Q* - I)) 2 

x [4(1 - 3u>(-i + n* A ) + n* A )(2 + 3w(n* A - 1) + 3n* A )+ 

+ 9e 4 n t 4 (l + w) 2 (n* A - 1)(2ST A + 1)(1 + w(2ft* - 1)) + 
+ 6e 2 t 2 (l + iu)((8 - 13Q A )n* A + 3u> 2 (ftl - l) 2 (2fi^ - 1) + 

-2w(n* A -i)(-2 + n* A (5 + n* A )))} 

r 2 

(1 + w)(Q* A - l) 2 {6w{Q* A - 1) - 2(Q* A + l) + 3e 2 t 2 (l + w)(2Q* A - l)) 2 

x [4(-i + 3w(n* A - 1) - + - 1) + ^a)(-i+ 

+3w(fi; - 1) + 3Q* A ) - 9e^ 4 (l + w) 2 (Q* A - 1)(2Q* A - 1)(2Q* A + 

+w(-i + w(n* A - i)(4n* A - 3) + 01(4^1 - 1)))+ 

+6e 2 t 2 (l + w){3w 3 (Q* A - If + w 2 {Q* A - 1) 2 {Q* A - 1)+ 

+(Q* A + 1)(3Q* A - l)(4ft* - 3) + iu(l + Q* A (9 + Q* A (27Q* A - 37))))] . 

The general solution of (50) can only be obtained by numerical analysis. Nevertheless, it 
admits an analytical solution in the form of a power-law <f>(t) = $ t p in the long wavelength 
limit e n — > 0. By direct substitution, one can find the two solutions 

2 4 

p - = 1 ~ (i + w)(i-n* A ) ' p+ = 2 ~ 3(i + w)(i-n* A ) (52) 

that do not depend on the values of A* and g*, but only on the physical quantity Q* A , as it 
was expected. 

While the value of the mode p_ is always negative (remember that in a spatially flat 
Friedmann- Robertson- Walker < Q* A < 1), p + can be either a growing or a decreasing 
mode. In particular, it changes sign from positive to negative for 

n* A = 1 - — — — - (53) 
A 3(1 + w) K ' 

that becomes Q* A = 1/3 in a matter-dominated regime and Q A = 1/2 in radiation-dominated 

one. This happens to be the same range for which the deceleration parameter q is negative 

(see section II). This means that in the present framework an accelerated expansion of the 

background, with q < and thereby a > 1, corresponds to a decreasing of the perturbations 

ip± < 0), as expected in a standard inflationary scenario. 
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The argument just described is a surprising similarity between the outcomes of the RG 
improved cosmology and those of the well-known inflationary theories based on the existence 
of an inflaton field [34-36], characterized by a non- vanishing Vacuum Expectation Value that 
plays the role of a vacuum energy. Adding this up to what already presented in [18], one 
can conclude that up to now RG improved cosmology is a self-consistent alternative to the 
usual inflationary paradigm. 

V. QUASI-CLASSICAL EVOLUTION OF PERTURBATIONS 

It is now worth checking if the RG improved cosmological perturbations show a late time 
behaviour that is in agreement with the usual ACDM scenario. In order to achieve this, one 
has to apply the evolution equation (48) to the background described at the end of section 
II. The starting point is the trajectory descibed in equation (15), that can be put in the 
closed form 




(54) 



which can be inverted to give the relation 




(55) 



where e = ±1 in the branch of the trajectory where k ^ kx- 

For what regards the Hubble rate, its evolution equation with respect to A is 




(56) 



and for the /3-function one can write 




(57) 



so that (56) admits the solution 




(58) 



where H T = H(\ T ). 

The last quantity that needs to be computed is the matter energy density 




(59) 
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that using (55) and (58) becomes 



32tt 2 " V 3 ; A + ev ^^ 



2 



One can simplify the expressions above obtained using the variable z = A/At, getting for 
the dimensionless Newton constant 

g(z) = 4ttA t (z + eV* 2 - l) (61) 

for the Hubble rate 

if (*) =ff r U + - 1J (62) 
and for the matter energy density 

3H T c2 ( e*T \ (z + v^rf* 

Plugging the above results into equation (48) one ends up with the following second order 
differential equation for the gradient 

AD li + BD IM + CD li = (64) 

where A, B and C are z-dependent coefficients. In order to solve the (64) in terms of the 

new variable z, one must keep in mind that, given a function / 

df_ _ d^chdXdk _d£ 1 dk _df 1 1 dH _ 
dt dz dX dk dt dz At k dt dz At H dt 



2 V ' \ 3 J V 2ttA t / cte 
so that, being the harmonic expansion of = ra \I/™, equation (64) becomes 

C 2 (z) <j>»(z) + C 1 (z) f\z) + C (z) <f>(z) = (66) 

where the tilde was used to underline the fact that, due to (65), the coefficients have changed. 

Equation (66) admits for z ^> 1 an analytical solution in the form of a power law. In this 
limit, the coefficients assume the form 

9{l + wfQ T 



1 + 4e + 3w(l + 2e) 

_L A A 



C (z) ~ 9( 1 + ^) 3 ^6(3 + 46) 

C1{Z) - (l + e)(l + 4e + 3w(l + 2e)) * (6?) 

r t \ 18(1 + W ) 3 ^t6 2 3 , Ra x 

Co( ^ (i +e )(i + 4 f+ 3m + 2 £)) z (68) 
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where £It = £ 2 At/3. Writing = Az a , one finds for the exponent the (w-independent) 
values 

a = and a = (69) 

that for the positive branch (i.e. k > kr) means = —2 and a + = —1/2, while for the 
negative one means «_ = and a + = 1/2. For what regards the positive branch, one has 
to make clear that the dimensionless cosmological constant is decreasing with respect to 
the cosmological time, so the perturbations increase. The negative branch results are even 
more interesting, because they describe the evolution of the perturbations during the quasi- 
classical regime of the RG improved cosmology. In the limit for k <C kx the cosmological 
constant can be neglected and the background can be approximated to a standard Einstein- 
Hilbert (EH) cosmology, so the evolution of its scale factor can be described to a good 
accuracy as 

a(t) = a t^^ (70) 

so that inverting the relation gives 

3(1+™) 



t(a)=[±) 2 (71) 



and beacuse 



one ends up with the result 



"M-aoW < 72) 



2 



H(a) = — — (73) 

V ; 3(1 + w) \aoJ 



but from (15) it is known that for k <^ k T 

A 

X r '~2k 2 2 I! 



A \k\ 1H T 
^=^^ot| = oT7| (74) 



so that it can be written 

9(1 + ™)W« X3(1+W) 



z(a) = v - ' H T ( - ) (75) 



and then the growing mode of the density perturbation goes like 

3(l+tu)/2 

(n\ - A- I — } 

y grow 



a) = 0o ( - ) • (76) 

00/ 



For a radiation-dominated universe one has that the growing mode grO w goes like a 2 , 
while in a matter-dominated one it goes like a?/ 2 . This is only partially in agreement 
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with the standard theory proposed in [37]: from there the evolutions were found to be 
oca 2 and oca respectively, so that only the radiation-dominated regime coincides with the 
present analysis. Nonetheless, one has to keep in mind that the present realization of the 
renormalization group is based on an averaging procedure over momentum scales greater 
than H, so that it only makes sense to give a position dependence to the coupling constants 
when they are evaluated over distances £ > if -1 . All the modes lying inside the curvature 
radius should therefore be evaluated at the averaged cutoff scale H, so that the spatial 
gradients T M and A M defined in (33) should automatically vanish. This is true in particular 
for the very late stages of the universe, when the modes of the perturbations that lie inside 
the observable window reenter the horizon. Such statement is consistent with the long 
wavelenght approximation used throughout the paper, as for spatial momenta greater than 
the Hubble rate the formalism should break down and a switch to the standard evolution 
should be mandatory. 

VI. CONCLUSIONS 

The analysis performed, by means of the gauge-invariant cosmological perturbation the- 
ory, provided a description of the evolution of the energy density fluctuations around a 
background where the coupling constants vary according to renormalization group equa- 
tions. 

Analytical solutions were found only for perturbations with very long wavelength com- 
pared to the typical curvature radius of the backgrounf universe. An intriguing issue would 
be to study the wavelength dependence, in order to construct a power spectrum of the 
primordial cosmological perturbations. This will be the aim of a future communication. 

A special thank to Roberto Percacci, Carlo Baccigalupi and Alfio Bonanno for useful 
suggestions and comments. 
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